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Abstract 

We discuss non-Gaussianity and baryonic isocurvature fluctuations in the curva- 
ton scenario, assuming that the baryon asymmetry of the universe originates only 
from the decay products of the inflaton. When large non-Gaussianity is realized in 
such a scenario, non-vanishing baryonic isocurvature fluctuations can also be gener- 
ated unless the baryogenesis occurs after the decay of the curvaton. We calculate 
the non-linearity parameter /nl and the baryonic isocurvature fluctuations, taking 
account of the primordial fluctuations of both the inflaton and the curvaton. We 
show that, although current constraints on isocurvature fluctuations are severe, the 
non-linearity parameter can be large as /nl ~ 0(10 — 100) without conflicting with 
the constraints. 



1 Introduction 



Current cosmological observations are now very precise to give us much information about 
the early universe. In particular, from observations of the cosmic microwave background 
(CMB) such as WMAP [11 [2], we can probe the physics of the early universe through 
the nature of primordial fluctuations, which are usually characterized by the amplitude of 
the fluctuations and its scale dependence as well as by a possible contribution of gravity 
waves. Recently in addition to these quantities, non-Gaussianity has been attracting much 
attention since it provides information about different aspects of the physics of the early 
universe. Non-Gaussianity is usually quantified by the so-called non-linearity parameter 
/nl and the recent constraint on this quantity from WMAP5 is [ij^ 

-9</nl<111 (95%C.L). (1) 

Although the data is still consistent with Gaussian fluctuation, which corresponds to 
/nl = 0, the central value is away from zero as /nl ~ 50. If future observations confirm 
such non-zero large value of /nl, it will give very important implications to the scenario of 
the early universe; a simple model of inflation would be excluded since fluctuations from 
the inflaton are almost Gaussian. In such a case, some mechanism is needed to generate 
large non-Gaussian fluctuations. 

Importantly, the simple inflationary scenario is not the only possibility to generate 
the density fluctuations. In particular, from the viewpoint of particle physics, there may 
exist a scalar field other than the inflaton, i.e., so-called the curvaton [31 IH Ej, which 
acquires primordial fluctuation and generate the present density fluctuations although it 
is a sub-dominant component during inflation. With the curvaton, large non-Gaussianity 



can also be generated [6l [7] |^ and the curvaton scenario seems attractive in the light of 
constructing a successful model of generating large non-Gaussianity. 

In order to generate large non-Gaussianity with the curvaton, it is necessary that the 
energy density of the curvaton at the time of its decay should be much smaller than that 
of the dominant component of the universe (which is expected to be from the inflaton). 
Density fluctuations generated in such a scenario are (almost) scale-invariant, and hence 
can be consistent with the observations if the fluctuations of all the components are adia- 
batic. However, in such a scenario, there exist entropy fluctuations between components 
from the inflaton and those from the curvaton. In order to have vanishing isocurvature 
fluctuations, it is necessary to generate baryon asymmetry and cold dark matter (CDM) 
at low-temperature universe after the decay of the curvaton. Although there are several 
possibilities that the dark-matter density is determined at a relatively low temperature, 
such as axion dark matter, the lightest superparticle dark matter and so on, baryogenesis 



'^-'^Here /nl represents the so-called "local type" non-Gaussianity. In this letter, we only consider non- 
Gaussianity of this type. 

"^^Some other models generating large non-Gaussianity have also been discussed, such as the modulated 
reheating scenario [8|. 
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at a low temperature may be challenging. Indeed, for some of the scenarios of baryo- 
genesis, like the thermal P and non-thermal [10] leptogenesis, a relatively high cosmic 
temperature is required!^ This fact indicates that, to realize large non-Gaussianity, it 
may be necessary to generate the baryon asymmetry of the universe before the decay of the 
curvaton. Since the size of baryonic isocurvature fluctuations is now severely constrained 
by observations, such constraints are important in considering baryogenesis in the cur- 
vaton scenario generating large non-Gaussianity. In particular, in the simplest curvaton 
scenario where the cosmic density fluctuations are totally from the primordial fluctuation 
of the curvaton, large non-Gaussianity cannot be generated without conflicting with the 
constraints from the baryonic isocurvature fluctuations if the baryon asymmetry is solely 
from the decay product of the inflaton. This argument excludes some of the scenarios of 
baryogenesis which require high cosmic temperature in the curvaton scenario. 

However, in the mixed fluctuation scenario where fluctuations from the inflaton and 
the curvaton both contribute to cosmic density fluctuations [HIIIS], the situation changes; 
even if the baryon number of the universe originates only from the decay products of 
the inflaton, the amplitude of isocurvature fluctuations relative to adiabatic ones can 
be suppressed. Thus, in such a case, large non-Gaussianity may be generated without 
conflicting with the constraint on the isocurvature fluctuations. 

In this letter, we discuss non-Gaussianity and baryonic isocurvature fluctuations in the 
curvaton scenario where only the decay products of the inflaton are responsible for the 
baryogenesis. We pay particular attention to the question how large /nl can be while 
imposing the isocurvature constraints. We will see that /nl ~ C(10 — 100) can be realized 
even if severe constraints on isocurvature fluctuations are imposed. 

2 Scenario and Formalism 

We first describe the scenario we consider in this letter and summarize formulae to calculate 
density fluctuations and its non-Gaussianity. 

Here, we consider the curvaton scenario where fluctuations from the inflaton also con- 
tribute to the total curvature fluctuations along with that from the curvaton. In addition, 
we assume that only the decay products of the inflaton (not those of curvaton) are respon- 
sible for the baryon asymmetry of the universe; the possibilities include a baryogenesis in 
the thermal bath before the decay of the curvaton and that in association with the decay 
of the inflaton. We also assume that the density of CDM is determined after the decay of 
the curvaton so that there is no CDM isocurvature fluctuations. We can easily apply our 
results to the case where CDM originates only from the decay products of the inflaton. 

In our study, we adopt the following form for the scalar-field potential: 



T^^However, the electrowcak baryogenesis [TT] may be a possibility. Another possible model is AfHcck- 
Dine baryogenesis [H] and the issues of non-Gaussianity in the model is discussed in [13j . 




(2) 
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where and a represent the inflaton and curvaton fields, respectively, V{(j)) is the potential 
for the infiaton, and rrifj is the mass of the curvaton. (We assume a quadratic potential 
for the curvaton.) In general, the curvature fiuctuations from the infiaton depend on the 
infiaton potential. However, for the purpose of the following discussion, we only need to 
specify the value of the slow-roll parameter e, which is defined as 

l.,2 /V, 



,M'Af], (3) 



where = dV/dcj), and Mpi ~ 2.4 x 10^^ GeV is the reduced Planck scaleCJ We 
also assume that the mass of the curvaton is small and the curvaton acquires primordial 
fiuctuation during the infiation. (We denote the initial amplitude of the curvaton as a*; 
hereafter, the subscript "*" is for quantities at the time of the horizon exit.) 

After the infiation, the infiaton begins to oscillate around the minimum of its potential 
and then decays. We denote the decay rate of the infiaton as and define the reheating 
temperature Tr as 

/ 10 

Tr=[^,MItI] , (4) 

where we use ^^sm = 106.75 as the effective number of the massless degrees of freedom. As 
the universe expands, the expansion rate of the universe H becomes comparable to rrifj 
and the curvaton starts to oscillate. Then, when the expansion rate becomes comparable 
to the decay rate of a, which is denoted as Fq- (and is related to the lifetime of the curvaton 
as To- = F~^), the curvaton decays. The start of the oscillation and the decay may occur 
before or after the reheating due to the inflaton decay, depending on the values of mg. and 
F.. 

Assuming that the potential of the inflaton is well approximated by a quadratic one 
around its minimum, its energy density behaves as that of matter for the period of the 
inflaton oscillation. Then, denoting the energy densities of radiation components from the 
decays of (j) and a as and p^^, respectively, evolutions of these variables (as well as 
those of the energy density of the inflaton fleld and the curvaton amplitude) are governed 
by 

p^^ + ^Hp^^ = F^p^, (5) 
p^a + AHp^^ = F^cr^ (6) 
P0 + 3Hp^ = -F<^p0, (7) 
a+i3H + T„)a + mla = 0, (8) 

where the dot represents derivative with respect to the cosmic time and pi indicates the 
background energy density of the component i. Notice that, when H -C m^, Eqs. ([6]) and 



T^^In the following, we assume that the primordial fluctuations are (almost) scale-invariant. 
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([8]) can be well approximated by 



P7. + 4i/p7, = r<,p^, (9) 
p„ + 'iHp„ = -T„p„, (10) 

respectively, where pa is the energy density of the curvaton field. In addition, it should be 
noted that the total radiation energy density is given by 

Pi = Pi^+Pi.- (11) 

In our analysis, we use the 5N formalism [TB] to calculate the perturbations. Then, 
perturbation variables are obtained by evaluating the number of e-folds from some time 
during infiation to the time after the curvaton decay as a function of model parameters 
such as T fj^rrifj and a*. Here, let us briefiy summarize the resultant formulae of the 
density fiuctuations and non-linearity parameter /nl- 

With the 5N formalism, if there is no isocurvature fiuctuation, the curvature fiuctuation 
originating from scalar-field fiuctuations is given by 

C^'^'^') = NJ^l + ]^Nab5^:5^l + ■ ■ ■ , (12) 

where N is the number of e-folds, 5ip1 is the primordial fiuctuation of the scalar field v?'*, 
and 



dN d^N , , 



with ip°' = (f) and a in our case. 

In the following, we consider the scenario where fiuctuations of the infiaton and the 
curvaton are both responsible for cosmic density fiuctuations. For simplicity, we assume 
that these fields are uncorrelated. Then, the curvature perturbations originating from 
these scalar fields can be calculated separately. The curvature fiuctuation from the infiaton, 
which we denote C^'^'^''*; is given by 

cr^ ^ -7^^<p* + i fi - ^) ^<pi (14) 



v^Mpi 2 V 2e 

where rj = M^]V^^/V, and we have used the slow-roll approximation. (We neglect terms 
of the order of 6(f)l which are irrelevant for our discussion.) In addition, the curvaton 
contribution to C*-'^'^'^ which we denote is expressed as 

^(adi) ^ j^j^^ ^ ^N^Jal (15) 

We solve the set of equations Eqs. ([5]) - ([8]) with parameters F^, F^., nia and a^, being 
fixed, then we determine N^j and N^^a- to obtain For the important case where 
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\Pia / P'^]t^r„ ^ 1 (^^d rrifj <^ r^), in which large non-Gaussianity can be generated in the 
curvaton scenario, the relation 4A^(cr*) = [P7„/P7]t>T„ holds and the a* dependence of 
is well approximated as [T5] 



3V27r 



r2(5/4) 



for [p^iJ p^]f»r^ < 1, 



(16) 



while, for the case where the curvaton eventually dominates the universe, [P7^/P7]t>Ta — 1; 
it is given by 



for [piJp^]f»T^ - 1- 



(17) 



Now, we discuss the entropy fluctuation between baryon and radiation, which is given 
up to the second order by 



5pb 1 / 5p, 



Pb 2 V Pfe 
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6p^ 1 f 6p. 



Pt 2 V P7 



Here, 6pi denotes the energy- density fluctuation of the component i on the uniform density 
slicing. We evaluate Sb^y when '^t <^ icq with tgq being the time of the radiation-matter 
equality. In the present scenario, there exist two sources of radiation: one originating from 
the inflaton and the other from the curvaton. We treat them separately and write 5p^/ p^ 



as 



5p^ 

P7 



^Pis + ^P 



7<T 



P70 + Pl^ 



(19) 



Since the baryon asymmetry and 7^ has the same source (i.e., the inflaton), the adiabatic 
relation holds between these two components. Furthermore, as will be discussed later. 



large non-Gaussianity can be generated when [p^„/ p^ 



7J t>T<^ 



<^ 1, thus we concentrate on 



such a case. Neglecting the terms which are second or higher order in [P7<^/P7]i^^^, 
using the adiabatic relation between the baryon and 7^, we obtain 



and 



Si 



67 



Pt. 
Pi 



(20) 



where terms of the order of {S p^y^ / Pj^)"^ vanish. (Here and hereafter, it should be un- 
derstood that [P7^/P7]t>r<^ is equal to 4A^(cr*), and is proportional to a^.) Notice that, 
compared to 6p^^/p^^, (5p^^/p^^ is of the order of [P7^/P7]t-^^_^ and hence its contribution 
is irrelevant in the present discussion. 



When [p7,/p. 



<^ 1, the cosmic expansion is solely determined by radiation from 
the inflaton, and we can neglect the effect of 7o- on the background evolution. Then, p^^ 
is proportional to and 



which results in 



b-y 



' Pi.' 




.Pi . 





2 a? 



(22) 



In discussing the non-Gaussianity in the present framework, it should be noted that 
the isocurvature fluctuations also generate curvature perturbation. Indeed, using the re- 
lation SC*-^™-* = 5*^7 (with Sm-y being the entropy fluctuation between the total matter 
and radiation), which holds in the matter-dominated universe, we obtain the isocurvature 
contribution to the curvature fluctuation as 



Scr on 





'Pl.' 




2 


.Pi . 





6a ^ 



l_5al 
2 al 



(23) 



where VLb and VLm are density parameters of the total matter and baryon, respectively. In 
our numerical analysis, we use Vti^/Vtrn — 0.17 [Ij. 

Since the fluctuations of the inflaton and the curvaton are assumed to be uncorrelated, 
the amplitude of the total curvature fluctuation C, in matter dominated epoch is given by 



c = 

^^adi) _|_ ^^iso) 



(24) 



where C,^ = C^'^^^ and Co- = + (a""' ■ From the observations of the cosmic density 
fluctuations, the size of ( is constrained. In our study, we determine the amplitudes of the 
primordial scalar-fleld fluctuations so that ( becomes consistent with the observed value. 

Non-Gaussianity of fluctuations is usually quantifled with higher order statistics such 
as bispectrum. Here, we consider the bispectrum of (: 



Then, B(^{ki, k2, k^) is obtained as 



(25) 



Bdku h) =[K + KAL In k„,i,,L N,, [PsAki)P5a{k2) + (2 perms.)] , (26) 



where 



N.rr 



2 



2 ^^T-ri 



VPi J 

Pi. 
Pi 



Here, Px{k) denotes the power spectrum of the variable X deflned as 

{X^,X^,,^) = i27rf6ih + k2)Pxik,), 



(27) 
(28) 

(29) 
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and is related to as 



(30) 



Assuming that the curvaton fluctuation is due to the quantum fluctuation during inflation, 
we obtain 



(31) 



where is the expansion rate during inflation and we assume A^^ to be scale- invariant. 

In deriving Eq. ( 126|) , following [171 CSj , we have regularized the infrared divergence by 
introducing the the infrared cutoff parameter and fcmin = min(/ci, /c2, ^s)- (We have 
neglected an 0(1) coefficient in front of Infcmin-L.) In order to discuss the implication to 
the cosmological observations, both and are taken to be the cosmological scale. 
Since the scale dependence from In k^^L is rather weak, we approximate In k^i-aL = 1 in 
our numerical analysis. 

As mentioned in the introduction, the non-linearity parameter /nl is often used to 



characterize non-Gaussianity of fluctuations, which is given 



Bcik,, k2, k;) = pNL [Pc(^i)^c(^2) + (2 perms.)] . 
Here P(^{k) is the the power spectrum of ( defined as 

In the present scenario, P(;{k) becomes 



A^ 



27r2 



2eMl 



A 



Ar 

ua do 



In kL 



oo 



(32) 



(33) 



(34) 



Then, the non-linearity parameter /nl is given by 



6 , _ (26 - r^) + Ae'M'^,{Nl + N^^l In k^^L )N^ 
1 + 26M2 (iV2 + Nl^Al In k^^L) 



fi 



(35) 



Since there exists the contribution from isocurvature fluctuations in this scenario, /nl here is not 
the same as the one for the case only with adiabatic perturbations. However, as shown in Eq. (|23p . the 
isocurvature contribution is suppressed by the factor flb/^lm- In fact, because of the difference between 
the transfer function for the adiabatic contribution and that for the isocurvature one, the bispectrum is 
enhanced at the Sachs- Wolfe plateau. In the present case, the enhancement factor for the corresponding 
term is roughly estimated to be ~ 2^/'^ for large-scale fluctuations [19]. Even if the second terms in 
Eqs. ([27]) and ((28l) are multiplied by this factor, the contribution from the isocurvature fluctuations is still 
sub-dominant. For the details of non-Gaussianity from isocurvature fluctuations, see [l9] . 
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where we have used the relation A^^ = A^^. Since we are interested in the case where 
the non-hnearity parameter becomes large, we neglect (2e — 77) in the numerator of (135|1 
in the following analysis. We can find approximated formulae of No- and iVg-o- for the most 
important case of [p-y^ / p'y]t:^T„ -C 1, for which /nl ^ 1 may be realized. In such a case, 
as we have mentioned, the relation AN{a^) = [p-y^ / P-y]t^T^ holds and hence 

N^, N.. - (1 - • [PiJPi\t»r. « 1- (36) 

For the sake of the following arguments, we also calculate the power spectrum of the 
isocurvature perturbations. With Eq. fl22l) . we obtain 



Ai = - 

Obi 4 



'pl.' 


' { 


.Pi . 





^ + ^Al\nk^i^L]Al. (37) 



Before showing the numerical results, we briefly consider the pure curvaton case where 
all the cosmic density fluctuations are only from the curvaton. In such a case, we obtain 



As,, N, 



acr 



Ac K IP,} 

and As^^ and A^ are of the same order irrespective of a*. (See, for example, Eqs. ( |T6l) 
and (ITTl) .) In this case, the entropy fluctuation is too large to be consistent with the 
observations [5], El EQl El] . Thus when the baryon number is produced at high temperature, 
it is difficult to have large non-Gaussianity in the simplest curvaton paradigm without 
conflicting with the isocurvature constraint. As we will see in the following, the situation 
changes in the mixed fluctuation scenario. In particular, when the curvature perturbation 
mainly comes from the inflaton fluctuation and [p-y^ / P-yjt^r^ ^ 1, large non-Gaussianity 
becomes possible without conflicting with the isocurvature constraint. 



Pya 



(3^ 



3 Numerical Results 

Now, we show our numerical results. In our analysis, we numerically solve Eqs. (jSD ^ dH]) 
and calculate the number of e- folds as a function of a,,. Then, we calculate /nl and Siyy 
for various values of the model parameters. In the following, we take m„ = 100 GeV. 

First, we show how large the baryonic isocurvature fluctuations can be. In Fig. [1], 
we show contours of constant As^^/A(^ on the a^, vs. T^/m^ plane. Here, we take Tr = 



10^° GeV, and e = IO-^e!! When > N^^Aj^lnk^i^L, one can see that the baryonic 



isocurvature fluctuations are suppressed as a* becomes smaller. This fact can be easily 



"^^In our following numerical analysis, we use e = 10~^ and 10"^" for illustrational purposes. In fact, the 
value of e = 10~^ corresponds to the case of quadratic chaotic inflation with Nf. = 50, while e ~ 0(10^^°) 
is realized in some class of new inflation model. 
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]^g-l4 r I il I I I I I null 

10-8 lO-'' 10-6 lQ-4 10-3 10-2 

CT. / Mp, 

Figure 1: Contours of constant A5j^^/A^ on vs. Ta/ma plane for Tr = 10^° GeV. For the 
slow-roll parameter, we take € = 10"2. Re gions excluded by the constraints (f40]l are shaded. 



understood as follows. When a* is small enough, the curvature perturbation is dominated 
by the inflaton contribution (,4, and hence is independent of a*. In addition, in such a 
case, the ratio [p-^/^ / P'y]t^Ta is proportional to a1 and hence Ag^^^/A^ becomes smaller as 
(T* decreases as far as the first term in the parenthesis of Eq. ( 1371) dominates. On the 
contrary, for N^. < iVj^ A^^ In fcmin-^, the isocurvature perturbation A^^^ is determined by 
the second term in Eq. (1371) . Then, A^^^^/A^ becomes insensitive to cr,., as shown in the 
figure. 

From current cosmological observations, baryonic isocurvature fluctuations are severely 
constrained; there is no sign of the isocurvature fluctuations in the observed angular power 
spectrum of the CMB, and Sh-y is consistent with zero. In our analysis, we adopt the bounds 
on the baryonic isocurvature fluctuations obtained from the latest WMAP5 result. We 
classify the baryonic isocurvature fluctuations into correlated and uncorrelated parts as 



[5, 



- As,^ sin 5, 



[5, 



A5, 



67 



1 — sin 5, 



(39) 



where sin 5 = A^^/A^. Then, we adopt the following bounds on the ratios of these 
isocurvature modes to A^, reading off the 95 % C.L. constraints from the WMAP5 results 



Si 



67 



A 



CJ 



> -0.31, 



S, 



A, 



< 1.35. 



(40) 



Notice that, in the present scenario, [S'f,^/A^]corr is negative^ In Fig. [U we shaded the 



'^'^The bounds given in are separately obtained for the cases of totally correlated and totally 
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10-8 lO-'' 10-6 10-5 10-* 10-3 10-2 10-12 10-11 10-10 10-9 10-B 10-7 IQ-e 

a, / Mp, CT, / Mp, 



Figure 2: Contours of constant /nl on the a* vs. Ta/ma plane for Tr = 10^° GeV. (/nl = 1, 
10, and 100 from the top.) For the slow-roh parameter, we take e = 10~^ (left) and lO"^*' (right). 
Regions excluded by the constraints (j^Oj) are shaded. 



region which is excluded by the above constraints; we found that the constraint on the 
correlated one is more stringent. 

It should be noted that, in the allowed region, the curvature perturbation ( is mainly 
from the inflaton fluctuation. Then, A^^ and A^^ are related to A^, which is known from 
the observation of the cosmic density fluctuations, as 

AL = A^^ = 2eM2 A^. (41) 

Neglecting the scale-dependence, we adopt A^ = 2.457 x 10^^ pLj to evaluate /nl given in 
Eq. (135]). 

Now we discuss how large /nl can be in the parameter region consistent with the 
constraints fHOl) . In Fig. [21 taking Tr = 10^° GeV, we show the contours of constant /nl 
for e = 10~^ and 10~^°. With such a choice of the reheating temperature, the curvaton field 
starts to oscillate after the reheating by the inflaton. If the curvaton begins to oscillate 
after the inflaton decay, the number of e-folds after the inflation depends only on the 
combination of T^j/m^ once a^, is fixed. Thus, we show our results in the a* vs. T^j/rrifj 
plane. (Notice that, for different values of Tr and nia, the figure is almost unchanged as 
far as -C F,^.) On the same figure, we shaded the region excluded by the constraints 

uncorrelated isocurvature fluctuations. Since the bounds for the case where the correlated and uncorrelated 
isocurvature fluctuations coexist are not available, we adopt the constraint (|40p as reference values. In 
addition, since the constraint on the correlated isocurvature fluctuations are not shown for [5'b^/A^]corr < 
in [1], we adopt the constraint for the case of [Sb^y / Aq]coyi > to derive the constraint (|40|) . assuming that 
the bound on | [S'b^/A^Jcon l does not signiflcantly depend on the sign. For the validity of this assumption, 
see, for example, [22] . 
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10-19 ^ 




10-35 ^ 




10-20 ^ 



10-36 ^ 



10-21 ^ 



Tn = lOSGeV, e= 10-= 1 



10-3-' i- 



T„ = lO^GeV, €= 10-1" I 



10-22 



10-38 



10-e IQ--' 10-6 10-6 10-* 10-3 10-2 

(7. / Mp, 



10-12 10-11 10-10 10-9 10-8 10-7 10-6 
CT. / Mp, 



Figure 3: Same as Fig. [21 except for Tr = 10^ GeV. 



f l40|) . We can see that, even after imposing the isocurvature constraint, a very large value 
of /nl (i-e., /nl ~ C(10 — 100)) is possible with small enough a^; in such a case, even 
though the components related to the curvaton are always sub-dominant, /nl becomes 
large. (A possibility where a sub-dominant component generates large non-Gaussianity 
was first considered in [IE]-) In addition, when a* becomes small enough, /nl becomes 
insensitive to a*. This is because, in such a parameter region, the second term in the first 
parenthesis of Eq. (126!) . which is independent of a*, dominates. 

Next, let us consider the case where the curvaton decays before the reheating (i.e., 
during the infiaton oscillation). Even in such a case, large value of /nl can be obtained 
while satisfying the isocurvature constraints. To see this, in Fig. [31 we show the contours 
of constant /nl on cr* vs. T^/m^ plane for the case of = 10^ GeV. As one can see, in 
this case, /nl ~ 0(10 — 100) can be realized even if we impose the baryonic isocurvature 
constraints. 

4 Implications for Scenarios of Baryogenesis 

So far, we have seen that a large value of the non-linearity parameter of /nl ~ O{10 — 100) 
can be realized without conflicting with the baryonic isocurvature constraints, assuming 
Tji, a^, and as free parameters. If we fix the scenario of baryogenesis, however, it 
is often the case that a lower bound on Tr is imposed to generate large enough baryon 
asymmetry of the universe. So, finally, we discuss whether large /nl is possible for several 
scenarios of baryogenesis. 

In order to realize the thermal leptogenesis j9], Tf> should be higher than 10^~^° GeV 
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t5 Then, for = 100 GeV and e = lO'^, a* < 10^^ GeV and < 10"^ GeV are 
required to realize /nl = 10 without conflicting with the constraints (140!) . Notice that, 
even though To- has to be much smaller than mo-, the curvaton decay occurs when the 
cosmic temperature is T > 10^ GeV for = 10~^ GeV, which is well before the start of 
the big-bang nucleosynthesis (BBN)@ If rria is larger, /nl ~ 0(10) can be realized with 
a larger value of T^. 

Even though the scenario of the thermal leptogenesis is attractive, it is incompatible 
with some class of supersymmetric model because, if Tr > 10^~^° GeV, overproduction 
of the gravitino may occur [25]. If a scenario of baryogenesis which works at a lower 
temperature is needed, one of the possibilities is the non-thermal leptogenesis in which 
right-handed neutrinos are directly produced by the decay of a scalar condensation (like the 
infiaton) [10] . In the non-thermal leptogenesis scenario, the lower bound on is reduced, 
and is given by Tr > 10^ GeV [26]. As one can see, even with Tr ~ 10^ GeV, /nl can be 
as large as (9(10) (or larger) satisfying the baryonic isocurvature constraints. In this case, 
with = 100 GeV and e = lO'^, /nl = 10 requires a, < 10^^ GeV and < 10'^^ GeV, 
which corresponds to the decay temperature of the curvaton of ~ (9(10 GeV). 

In summary, even if the baryon asymmetry originates only from the decay products 
of the infiaton, a large non-Gaussianity of /nl ~ (9(10 — 100) is possible in large class of 
scenarios of baryogenesis without conflicting with the observational constraints. In such 
a scenario, however, baryonic isocurvature fluctuations are inevitably generated and may 
be just below the observational bound. Thus, if a large value of /nl is confirmed in future 
observations, it is strongly encouraged to look for signals of the isocurvature fluctuations 
to test baryogenesis models in the curvaton scenario. 

Note added: While we were finalizing the manuscript, Ref. [27] appeared on the arXiv, 
which has some overlap with our analysis. 
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